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1.  Introduction 

An  inqportwt  problem  in  understanding  the  properties  of  the 
nuclear  surface  is  the  natxire  of  the  deviations  from  the  Thomas-Fennl 
statistical  approximation  (TFA).  According  to  TFA^ the  particle  and 
energy  density  at  each  point  is  .-.hat  it  would  be  if  the  density  were 
uniform.  This  approximation  is  accurate  if  the  density  chfuiges 
sufficiently  slowly.  However,  in  the  region  of  the  nuclear  surface 
the  density  falls  from  (saiO  9DSS  to  10^  of  its  maximum  vsdue  in  a 
distance  conqparablc  with  other  relevant  distances,  i.e.  the  range  of 
nuclear  forces  and  the  Fermi  wavelength  of  the  system.  Thus  the  effect 
of  inhomogeneity  might  well  be  significant  here.  V/e  assxune  that  the 
energy  density  can  be  expanded  in  a  series  of  powers  of  <ll)  and 
powers  of  operating  on  S  ,  This  is  probably  a  serai-convergent 

expansion.  The  odd-  terms  vanish  by  reflectional  invariance  of  the 
Hamiltonian.  The  leading  or^er  terms  are  of  the  form  i{S)\P  S  \  ^ 

9 

♦  £(/)F  /•  The  latter  term  can,  after  partial  integration,  be  cast 

in  the  form  of  the  first  term;  thus  only  the  first  will  be  considered. 

A  correction  of  this  form  for  the  kinetic  energy  density  was  first 

proposed  by  WelzsAcker^\  However,  some  of  his  approximations  are 

2\ 

in  ejrror  and  more  recent  studies  '  indicate  that  the  inhomogeneity 
corrections  are  much  smaller  than  h'eizsacker  gave,  although  his  fora 
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is  preserved. 

In  this  paper  we  discuss  a  single  model  for  generating  density 
variations  for  v/hich  the  leading  order  inhomogeneity  effects  can  be 
calculated  explicity.  Instead  of  trying  to  study  the  behavior  of 
particles  at  the  nuclear  surface  explicitly  we  consider  a  set  of 
particles  moving  independently  in  a  cubical  box  upon  v*iich  is  super¬ 
imposed  a  sinusoidal  "ripple"  potential, 

^2 

u  “  ^ 

Although  this  is  a  different  situation  than  that  which 

exists  at  the  nuclear  surface,  it  is  possible  to  identify  uniquely 

2 

the  inhomogeneity  correction  terms  to  order  'if  p.  The  ripple 
model  has  aJ.so  been  considered  by  Lindliard''  in  his  study  of  the 
electron  gas. 


I 
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11.  CALCULAriOK  QF  rAliTICLL 

In  the  absence  of  the  ripple  potential,  tne  ave  function 
of  each  particle  is  .i  pLane  ^ave: 


ex? 


oxp  f  i  k^x  +  i  k  +  i  k^zj^(2) 


(  .’e  assume  periodic  boundary'’  condlt  one  and  a  box  of  unit  volume). 

The  presence  of  the  ripple  potential  chantcs  the  v;avefijnctio;-.  of  each 
particle;  ho'-'ever  the  total  wavef  ixct.  on  is  still  an  antis-TTimetrizod 
pro  luct  of  cinj le  particle  naval  lactions,  as  we  nctlcct  correlations 
due  to  two  body  interactions,  -ince  the  ripple  .lepends  only  on  the 
x-coorlinate,  the  y  and  z-dcp;.ndc  t  factors  ar-^.  unchani^ed.  ;*o’.;ever 
t'-.e  x-depender  t  art  of  cad;  nc’.  single  particle  ■.avefurxtion  is  a 
Ilathieu  function.  ..'e  can,  instead ,  express  the  new  wavef'anctions  in 
terms  of  the  ori^insQ.  plane  avos  by  a  perturbation  expansion.  Up  to 
second  order  in  the  potential  the  wavef  uriction  corresponding  to  an 
;n,..ert  urbed  ..-mentar.  k  is  ,  iven  as  folio;  s: 


(3) 


where 


-  6  - 


N  -  1  + 


fl  ,  1.  ■ 

^  Id?  d?. 


1  -il* 


and 


D^Oc)  -  nq(21^  +  nq), 


The  particle  density  in  state  k  is  Liven  by: 

2 

j>kW  -  *>1^-  1  -  ‘  CM  <?c  ♦  f-  ♦  3;^ 


-  e^3 


(4) 


To  calculate  the  total  particle  density,  vie  must  svun  over  aU 
occupied  single  pairticle  states  V/e  assume  that  the  particles  fill 
up  unperturbed  levels  In  order  of  increasing  energy  (neglect  any 
change  in  level  order  due  to  the  ripple).  Then  the  number  of 
particles,  neglecting  spin  amd  isotopic  spin,  is 

S  1,1  -  A  -  ^0  "  ^  '  (5) 

We  also  define: 

i>n  : 


1 
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;  ince  the  level  Icneity  i::  ej.a.ietric  1th  i-orji-ect  to  the  intcrchr.i.^^e 
follov/s  that  evaluation  of  the  term 

can  be  cieplified  by  use  of  the  iderhlties; 


“-1^-2 


n  [J - il 

‘■2=1  “a  ' 

1  r  1 

■  i?  ‘ 


rivar:  v’e  oVbain 


(7a) 

(7b) 


'£v^j^(x)  -  Ca^^coscx+  a^q'^r^  ^  cos2qx^(6) 

^  V\ 


;e:;t  we  calculate 


die 


(9) 


The  volame  element  is  -  inl^dl^dkj^ 

with  the  refion  of  inte£;:ation  for  Loing  from  o  to  (ky  -  k^)^,  and 


?  indicates  that  singularities  are  to  be  treated  in  the  sense  of 
principal  value, 

Jhus  e  obtain; 
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nqCek^  +  nq) 


(10) 


Now  introduce  new  variables: 


i’hen 


k  -  kp  X  j 

nq  -  akpy. 


P 


j 


♦  y) 


(11) 


A  straigthforward  calculation  gives: 


§ 


n 


^  \  +  ^fvT  1 1^1  n  ""tTo^ 

•-  ^  » ■»?  kk/ 


(12) 

IV 


In  the  limit  q  — ^o,  corresponding  to  a  long  wavelength  ripple,  the 
expression  for  the  density  becomes 


cosqx  ♦ 


cos2qx)« 


(13) 


This  agrees,  to  order  a  ,  with  the  Thonas-Fennl  result 


j^-pp(x)  “  (1  -  )(1  - 


cosqx) 


3/2 


(U) 
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ni.,  KINETIC  ENERGY 

V/e  now  turn  to  the  inhoniot®neity  correction  to  the 
kinetic  ener^*  The  expectation  value  of  the  kinetic  ener^ 
of  a  particle  in  state  k  is  given,  in  units  of  by 


- -K 

^  J  /w 


Substitution  of  (3)  into  (15)  yields 

2 


^  /  1  1  X 

K  ^D]Tk)  * 


T,.  -  k""  ♦ 


and  the  total  kinetic  energy  becomes 


(15) 


(16) 


(17.  a) 

-  2  2  2  -1 

i  -ft-  .  ft-a— t  A. 

(17.  b) 

q.^0  - 

This  result  agrees  to  order  a*  and  in  the  limit  with 
the  Thomas-Fermi  approximation,  according  to  v/nich 

^IF  -  V  K  (1 

-  5^^  *  D  *  f  <  ^>  ‘  f  <  (J?)*'  . 3  .  (18) 

where^^  denotes  average  over  space.  According  to  (13), 


cosqx  ♦ 


(19) 


.  --2*, 


v 


rhu! 


<^S>  vanishes,  but  ^  (6^)^ is  finite,  Tliis  gives 

T.jp  -  |kp  A(1  ♦  I  ^  ■*•  . . . ) , 


(20) 


in  agreement  with  (17b), 

In  a  non-unique  way,  it  is  possible  to  discuss  the  kinetic 
ener^  density, defined  by  the  requirement  Tj^  -  ^dr.  In  addition 
to  the  forms  suggested  £q. (15),  one  which  is  also  useful  in  the 
calculation  of  the  exchange  potential  energy  is 


Substitution  of  the  expression  for  the  wavefunction  (3)  into 
this  definition  leads  to 


^  2  ^  ^^x  *  ■‘'x  n^x 

t:^(x)  -  4  -  f-  (  - 5 - -  - ^)  - 


'J- 


'-1 


2^ 

+  —  |a  cosqx  ♦j^ 


fd^)^  ^ 


°1°2 


°-l®.2 


‘'x  1  a^  2 

- J  1“  C052(pc  ♦  (id 


"1  "-I 
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By  means  of  some  elementary  algebraic  manipulation  we  can  simplify  this 
formula  to 

0^3  fl  ♦  I  cos2qx;J 

(23) 

*  1  *  |cos2q3([J 

...)(1  +  icos2px)|^(24,) 

kp 

where  here  ^  y  denotes  average  over  the  distribution.  The  total 
kinetic  energy  is  the  spatial  integral  of  (24),  which  gives  Lq, (I7).* 

*  Although  the  average  kinetic  energy  density  is  given  correctly 

by  TFA^  the  explicit  "pression  is  not  quite  correct.  According  to  TFA, 
the  kinetic  energy  density  should  be  isotropic  at  each  point.  On  the 
other  hand,  according  to  our  ripple  model,  the  y- component  of  the 
kinetic  energy  density  is  i  iven  by 

while  the  x  component  is  not  simply  proportional  to  the  density.  This 
discrepancy  may  be  connected  with  >.he  fact  that  we  have  not  taken  into 
account  rearrangeiaent  of  the  level  order  due  to  the  ripple  potential; 
i.e.  the  levels  are  assumed  to  fill  up  in  order  of  increasing  unperturbed 
energy,  which  may  not  necessarily  ingly  order  of  increasinf  pepl^iy^y^  enerfy. 


k 


-  ^  cosqx  ^ 


/ar:riin£:  over  states  k,  vie  obtain 


k  ;  -  Tcosqx 

Z  & 


12 


In  order  to  interpret  the  nature  of  the  teims,  we 
reconsider  the  expression  for  6^ : 


(1 - ^  )  cos  (qx). 

121c: 


(25) 


•kj’ 


I 


Clearly 

s- 


2a  .  /  \ 

>mX 


(26) 


It  CA;:  c  i'cadllj'  verified  .t  the  chan^.e  in  hinetic  energy  due  to 
l.e  ripqlc  ccuils 

(2/) 


6T-  ^ 


The  first  tera  is  just  the  correction  discussed  above,  while  the 
second  term  is  an  inhoinoceneity  term.  This  correction  can  be 
expressed  as 


i  \SL£i 


*  So  (28) 


with 


I 


1 

9 


# 


The  v/cizsicker  correction 


1) 


of  just  tnis  form,  but  vdth 
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^  -  1,  Thus  our  correction  is  only  l/9  as  large  as  that  calculated  by 

V/eizsftcker 


It  ii,  ailEo  of  interest  to  consider  c-he  other  limit  of  a 
rat'i^ly  varying  potential,  i.e.  * 


In  this  case  we  have 


and  thus 


But  in  this  case 


so  that 


2a  cosqx 

- 2“^ 

q 


2a  sinqx 

q 


(29) 


(30) 


(31) 


(32) 


Thus  v.’e  find 

which  is  the  saoe  as  t'e  ’eitsachor  correction  “  1)«  However, 
V/eizs&cker '  s  original  derivat  on  was  intended  to  te  vsilld  in  the 
limit  q  — ^Oe 
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IV.  IMHO:-DGEMUirY  CQ.ufflCTION..:  TO  TI-Ii;  POrE-illAL  EIJlfiGY 

we  c6G.culate 

In  this  section /inhoino^eneity  corrections  to  the  potential 
energy  due  to  a  short  ran'^e  two  body  interaction,  but  neglect  all 
correlations  due  to  t’lese  interactions. 

First,  consider  a  local  (V/ij^ner)  interaction  and  neglect 
cxchanje  terms. 

The  total  interaction  energy  is  given  by 


(34) 


’i^nich  equals 


ijjf  • 


(35) 


Introducing  center-of-mass  and  relative  coordinates  U,  r,  we  can  write 
this  as 


(36) 


ilalcing  a  Taylor  exparision  of  the  density  about  its  value  at  R,  we 
find 

(a .  f)  -ym  .1  z  *••• 


(37) 
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To  calculate  the  potentiaQ.  energy,  we  assume  that  V  if  spherically 
symmetric.  Then  with  help  of  an  integration  by  parts,  we  obtain 

tS  I"  (38) 

The  first  term  represents  the  potential  energy  due  to  a  6-function 
interaction.  The  second  term  is  t  .e  lo. cst  order  classical 
inhorao^^  eneity  term  due  to  the  finite  range  of  the  nuclear  forces, 

.c;:t  consider  the  potential  energy  due  to  a  pure  space 
exchange  ^Majorana)  interaction.  This  gives 


ti.e  mixed  density  defined  by 

J)(ri,r2)  -  S 


(39) 


wnere 


(40) 


Again,  in  terms  of  a  center  of  -nass  and  relative  coordinates  we  obtain 


-  # 


■(R  +  §,  R  -  f )U(r)dhdr  . 


(41) 


Now  a  simple  expansion  about  H  shows  that 
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(42) 


whci-. 


-r^^(x)--X'  4*  -2  ^^-  ■  2  ^  ;i(43) 

T  ^  ^  5  X X ^  ^x”  0  Xj 

The  quantity  i®  closely  related  to  the  kinetic  energy  density 

defined  in  Eq.(24);  in  fact, 


X(3t)  -  ^i^'C^^Cx). 


(u) 


Thus,  for  a  centraQ.  potential  we  obtain 


-  I  ^dB  U<ij^  -  dfij  Ur^dr  ♦  9  (j^Ur^dr).  (45) 

The  volume  Integral  of^  can  be  calculated  using  (8)  and  (24).  It 
2 

is  seen  that  the  q  correction  2fO  times  as  large  as  the 

correction  to^X^*  Thus,  the  inhomogeneity  correction  to  ^j>T> 
ecMls  ^  \P_f  r  ’,  and  the  coirection  to  the  exchange  potential  energy 
is  only  1/9  as  large  as  for  a  direct  potential  of  the  same  strength 
and  range  ♦, 


*  An  exchange  Interaction  appears  In  this  respect  like  a  velocity  dependent 
6 -function  Interaction,  for  which  there  is  no  Inhoaogenel^  correction. 
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This  conclusion  appears  plausible  when  e  realize  that  the  -.matrix 
elements  of  an  exc.'.anf-e  interaction  depend  on  the  avora;  e  relative 
moment  '.. I  of  tlx  colliding  pair  before  and  after  the  scattering  but, 

(like  a  6 -function),  are  independent  of  the  momentum  transfer  in  the 
coldluion.  By  contrast,  the  matrix  elements  of  a  finite  range  local 
interaction  depend  only  on  the  momentum  transfer,  but  not  on  the 
momenta  t'.x..'iselves 
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V..  ONE  ai;d  r/.'o  diicncional  feri-h  gas 

It  is  possible  to  obtain  additional  insight  into  the  nature 
of  inhomogeneity  terms  by  considering  the  one  and  two  dimensional 
analogues  of  the  model  considered  above.  The  expressions  (8)  and 
(17. a)  still  hold  for  the  particle  and  kinetic  energy  densities. 
Hc'cver,  the  functions^  aa*e  now  ^  iven  by  different  expressions, 
since  the  integral  over  the  available  phase  space  is  different. 


In  the  one  dimen 

sional  case. 

we  have 

dk^ 

1 

“  ^ 

'  -kp, 

-kp 

dk.. 


(46) 


Making  the  same  change  of  variables  as  in  the  three  dimensional  case  , 
we  find 


i 


1 


(47) 


An  elementary  calculation  shows  that  in  this  case  the  inhomogeneity 

sign 

term  has  the  opposite  /ks  the  V.'eizsacker  correction  and  is  1/3  as 
large  i.e.  it  is  given  by  (28)  -ith^  **  “  3*  if  q^^Ekp, 

then  we  get  exactly  the  h'eizsacker  term,  th.e  same  as  in  the  three 


dimensional  case 


19 


In  the  two  dimensioned  case,  we  have 


It  is  :l^nlflclant  that  in  the  two  diner  ense,  we  obtain  no 

Inhomo^eneity  correction  o)  for  q  ^  2kp 

In  particidar,  wc  find,  for  q  ^  2kp, 


j)  -/oD  -;| 


with  no  term  of  order  a  , 


and 


(49) 


I 
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f  -5.  [I  <4  - 


2 

41^ 


♦  ^  (1  ♦  5  CM 


2^1 


(50) 


The  total  kinetic  enercy  is 

T  -  A  ^  ^), 

4k^ 


(51) 


-’’■dch  agrees  with  the  exact  TFA  value. 


It  is  also  interesting  to  note  that  in  the  t»«o  dimensional 
case,  the  density  is  proportional,  to  the  naxiinum  kinetic  energy.  This 
linearily  is  probably  closely  connected  vdth  the  absence  of  the 
inhomogeneity  term. 

\ 

We  vjould  now  like  to  investigate  i^aether  there  are  any 
correctiona  to  the  average  kinetic  energy  Involving  higher  pomre 
of  the  strength  g,  of  the  ripple  potential.  To  answer  this  question, 
we  must  calculate  higher  order  corrections  to  our  single  particle 
wavefunctlons.  It  can  be  demonstrated  expllclty  that  there  aure 

L 

no  further  corrections  to  the  kinetic  energy  through  order  a  , 
provided  I  iHit  we  have  not  verified  that  Eq.(5l)  is  exact  for 
two  dioensiona  (q^  2kp).  It  appears  (although  not  proven)  that 
corrections  of  order  a^  vanish  identically  provided 
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VI.  I.-UCLEAR  SURFACE  ^ 

If  the  form  of  the  inhomogeneity  correction  le  aesumeci  to  be 
of  the  form  (26)  in  nuclei^  then  an  analysis^  ^  shows  that  the  empirical 
▼alue  o£^  ,  including  both  kinetic  and  potential  effects,  lies  in  the 
range  o.36  to  o.72,  depending  upon  which  semi-empirlal  mass  formula  is 
used.  Since  .the  kinetic  contribution  is  only  1/9,  the  remednder  must 
be  attributed  to  potential  energy.  The  results  of  Section  IV  show 
that  thL  contribution  to  the  inho.mogeneity  of  an  exchange  potential 
is  only  1/9  that  of  an  identical  V/igner  potential.  Since  the  V/lgner 
and  I^Jorana  components  of  the  Interr.ucleonic '  potential  appear  to  be 
conparable,  the  major  part  of  the  Inhomogeneity  appears  to  be  derivable 
from  the  Vvlgner  psirt  of  the  potential,  l.e.,  it  is  a  classical  rather 
than  a  specifically  quantum  mechanical  effect. 

While  inhomogeneity  term  is  responsible  for  the  finite 
surface,  it  contributes  Just  half  (exactly  half  if  form  (28)  is 
assumed)  of  tr.e  surface  energy,  the  remaining  being  derived  from  the 
loss  of  binding  ener^  in  the  surface  region  as  would  be  computed 
from  (say)  a  Fermi  gas  model. 

Of  course  the  factor  1/9  is  valid  only  if  the  effective 
wavelength  of  the  density  variation  at  the  nuclear  surface  is  long 
compared  with  the  Fermi  wavelength,  that  is,  If  this 


bn«  applicability  of  a  ripple  model  to  the  nuclear  surface, 

estlinate  this  ratio  as  follows: 

Vfe  approximate  the  density  on  the  surface  region  by 
y  “  i  Jo  ^  ^  distance  0  Is 

related  to  q  by  qO  -  1,66,  thus  for  0  •  2,4  f  and  kp-  l,4f  ,  ws 
find  q/2kp  •  0.26,  vMch  appears  to  be  sufficiently  small  for  our 
purposes , 


In  conclusion,  we  reiterate  that  the  analysis  Is  mads  for 
an  independent  particle  model  vdth  local  (including  exchange) 
interactions .  A  more  comi^ete  analysis  oust  include  the  effect  of 
inter-nucleonic  correlations. 
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